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Reduced-order Distributed Consensus 
Controller Design via Edge Dynamics 

Dinh Hoa Nguyen, Member, IEEE 

Abstract —This paper proposes a novel approach to design 
fully distributed reduced-order consensus controllers for multi¬ 
agent systems (MASs) with identical general linear dynamics 
of agents. A new model namely edge dynamics representing 
the differences on connected agents’ states is first presented. 
Then the distributed consensus controller design is shown to be 
equivalent to the synthesis of a distributed stabilizing controller 
for this edge dynamics. Consequently, based on LQR approach, 
the globally optimal and locally optimal distributed stabilizing 
controller designs are proposed, of which the locally optimal 
distributed stabilizing design for the edge dynamics results in 
a fully distributed consensus controller for the MAS with no 
conservative bound on the coupling strength. This approach 
is then further developed to obtain reduced-order distributed 
consensus controllers for linear MASs. Finally, a numerical 
example is introduced to demonstrate the theoretical results. 

I. Introduction 

Multi-agent systems have gained much attention recently 
since there are a lot of practical applications, e.g., power grids, 
wireless sensor networks, transportation networks, systems 
biology, etc, can be formulated, analyzed, and synthesized 
under the framework of MASs. A key feature in MASs is 
the achievement of a global objective by performing local 
measurement and control at each agent and simultaneously 
collaborating among agents using that local information. 

One of the most important and intensively investigated 
issues in MASs (and their applications) is the consensus 
problem due to its attraction in both theoretical and applied 
aspects m, 121. The widely used consensus protocol ID for 
identical general linear agents has the following form 

Ui{t) = -fiK ^ = ( 1 ) 

where N is the number of agents, Ui{t) is the control input 
vector of Ah agent, Xi {t) is the state vector, jj, > 0 represents 
the coupling strength, K is a coupling matrix. A/) is the neigh¬ 
boring set of ith agent, > 0 are elements of the adjacency 
matrix (the definitions of graph-related terminologies will be 
given in Section III-Ab . 

Consequently, a natural question arose for the consensus 
law O is how to systematically design the coupling matrix K. 
One possibility is to employ the classical LQR optimal method 
of which K is the solution of an algebraic Riccati equation, 
e.g. a, a. Nevertheless, a only considered agents with 
integrator dynamics while 0 did not consider a global LQR 
performance index and hence the obtained distributed consen¬ 
sus controllers were locally optimal. The design of a globally 
optimal distributed consensus controller ([T]i for general multi¬ 
input multi-output (MIMO) linear MASs has recently been 
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presented in 16). Unfortunately, all the obtained results hitherto 
for the consensus control law O have conservative lower 
bounds for the coupling strength p. 

This paper firstly contributes a new approach to system¬ 
atically design distributed consensus controllers for linear 
MASs consisting of MIMO agents whose communications 
are represented by undirected graphs. The key idea is to 
derive a transformed representation of MASs that we call edge 
dynamics. This can be considered as a generalized notion and 
development of the introduced “edge Laplacian” in Q where 
only agents with single-integrator dynamics were treated. 
Based on the edge dynamics, we show the equivalence be¬ 
tween the distributed consensus controller design for the initial 
MAS and the distributed stabilizing controller design for the 
edge dynamics. Subsequently, two different methods namely 
globally optimal and locally optimal designs are proposed 
to obtain such distributed stabilizing controller. The globally 
optimal stabilizing design leads to a lower bound condition for 
the coupling strength of the consensus law, which coincides 
with other studies in literature. However, the interesting result 
for the locally optimal stabilizing synthesis reveals that the 
coupling strength can actually be arbitrarily small and positive 
such that the consensus is still achieved. To the best of our 
knowledge, this result on the consensus design has not been 
achieved so far. 

Next, the second contribution of this paper is on the design 
of reduced-order distributed consensus controllers. By further 
elaborating on the weighting matrices selection, the resulted 
LQR-based distributed consensus controller is shown to have a 
lower order than the previously obtained one. Especially, if all 
agents have a single pole at the origin, we show that the MAS 
can achieve consensus with a first-order distributed controller. 
Furthermore, all eigenvalues of the closed loop MAS and the 
consensus speed can be explicitly analytically determined. 

The organization of next parts in the paper is as follows. 
Section |II] introduces the edge dynamics and the design 
equivalence. Then systematic methods based on LQR approach 
to derive full-order and reduced-order distributed consensus 
controllers together with a numerical example are presented in 
Section|III]and Section ITVl respectively. Finally, some remarks 
are given in Section IV] 

11. Edge Dynamics and Design Equivalence 
A. Preliminaries 

Denote Q the graph representing the information structure in 
a multi-agent system, where each node in Q stands for an agent 
and each edge in Q represents the interconnection between two 
agents. In this paper, we assume that Q is undirected. Let us 
denote E £ incidence matrix and L the Laplacian 

matrix associated with Q, where M is the number of edges in 
the graph. Subsequently, we have L = EE"^. Denote the 
generalized pseudoinverse of the graph Laplacian L 0. The 
following lemma shows some useful properties of E and Lf 

Lemma 1: 0 Suppose that Q is connected then the fol¬ 

lowing statements hold. 

(i) = 0. 

(ii) LM = LLt = 
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Next, we introduce the notations used in this paper. The real, 
complex, complex with non-positive real part and complex 
with negative real part sets are denoted by R, C, C'L, and C_, 
respectively. The superscript * denotes the complex conjugate 
transpose. Next, 1„ represents the n x 1 vector with all 
elements equal to 1 , and /„ denotes the n x n identity 
matrix. Moreover, and ^ denote the positive definite and 
positive semidefinite of matrices, respectively. On the other 
hand, diag{} represents diagonal or block-diagonal matrices, 
Amin(^) denotes the eigenvalue with smallest, non-zero real 
part of A, and a (A) represents the eigenvalue set of A. Lastly, 
(g) stands for the Kronecker product. 

B. Edge Dynamics 

Consider an MAS consisting of N identical agents whose 
model is described by 

Xi{t) = Axi{t) + Bu^{t), ( 2 ) 

where A e B e Xi{t) G R" and u^{t) G R"" 

are the state and input vectors of the ith agent, respectively. 
Initially, the agents are independent, i.e., the whole MAS can 
be represented by 

x{t) = {In ® A)x{t) + {In ® B)u{t), (3) 

where x{t) = [xi{t)'^,... ,XN{t)'^]^, u{t) = 

..., UAr(f)^]^. Then each agent collaborates 
with some of other agents to achieve a specific global 
objective namely consensus defined as follows. 

Definition 1: The agents (|2]i are said to reach a consensus 
if the following condition is satisfied, 

lim \\xi{t) - Xj{t)\\ = 0\/ i,j = (4) 

i—>-oo 

Hereafter, the time index t is omitted in expressions of Xk, Uk, 
and other variables just for the conciseness of mathematical 
representations. 

The following assumptions are employed in this paper. 

Al: fT(A) G C^?., and at least one eigenvalue of A is on the 
imaginary axis. 

A2: {A, B) is controllable. 

The purpose of assumption Al is to avoid the consensus of 
agents to infinity or zero a whereas assumption A2 is for the 
existence of a controller. Let us denote a new state vector and 
a new control input vector 

Z = (g) In)x, W = {E'^ (g) Im)u, 

then multiplying both sides of Q with (g) In gives us 

i = (E'^ 0 In){lN 0 A)x + {E'^ 0 In){lN 0 B)u 

= (E^ 0 In)[{L^L + 0 A]x + (E^ 0 B)u 

= [{E'^L'^EE'^)®A]x+{Im®B)w. ( 5 ) 

Consequently, Eq. 0 can be further rewritten as follows, 

i = [(E^E^E) 0 A\z + {Im ® B)w. (6) 

We call the state-space representation 0 the “edge dynamics” 
of the MAS 0. The state vector z includes all the differences 


between the states of connected agents, i.e., it represents the 
MAS dynamics along the edges. Therefore, the initial MAS 
is consensus if this edge dynamics is stabilized. If each agent 
is an integrator, i.e., A = 0, E = 1, 0 becomes 


z = tu GG i = u, (7) 

which is much simpler than 0. The consensus protocol for 
this circumstance is u = —fiEx, which is equivalent to w = 
—gLLf.z, where Eg = E^E is called the edge Laplacian fT). 
Hence, the closed-loop dynamics of the integrator MAS is 


Z = —giLeZ -H- X = —fiLx. 


(8) 


The equivalent transformation 0 shows us a direct connection 
between the original integrator MAS and its transformed 
dynamics along the graph edges as well as the connection 
between the Laplacian and the edge Laplacian of a graph. 
More details can be found in El but are omitted here since 
they are unused in our consensus control designs. 

Denote E = E^E^E. The algebraic properties of E and Eg 
are presented in the following lemma. 

Lemma 2: The following statements hold if Q is connected. 

(i) Eg has exactly — 1 non-zero eigenvalues, which are 
equal to positive eigenvalues of E while all other eigen¬ 
values of Eg if exist are 0 . 

(ii) E has exactly N — 1 non-zero eigenvalues, which are all 
equal to 1, and other eigenvalues of E if exists are 0. 
Proof: (i): Let (A, v) be any eigen-pair of E with A 7 ^ 0. 

Consequently, AE^u = E^'^Lv = EgE^u, hence (A, E^u) 
is also an eigen-pair of Eg. If Q is connected then E has 
a single zero eigenvalue while all other eigenvalues of E are 
positive. Therefore, Eg has N—1 non-zero eigenvalues, which 
are equal to non-zero eigenvalues of E. On the other hand, 
trace(E) = trace(EE^) = trace(E^E) = trace(Ee). Since 
Eg ^ 0 , all other eigenvalues of Eg if exist must be 0 . 

(ii); Let (A, v) be any eigen-pair of E. We have XEv = 
ELv = LL^'^Ev = Ev. As a result, A = 1 is an eigenvalue of 
E. Any other eigenvalue of E if exists must have the associated 
eigenvector v satisfying Ev = 0 and hence Lv = 0 and EgU = 
0, i.e., that eigenvalue is 0. From the proof of part (i), the 
number of eigenvectors of Eg associated with 0 eigenvalues if 
exist is M — N + 1. Thus, E has exactly N—1 eigenvalues 
equal to 1 and all other eigenvalues if exist are 0 . ■ 

Let U G R^^^ be an orthogonal matrix that diagonalizes 
E, and z = {U^ 0 /n)z, w = {U'^ 0 Im)w. Subsequently, we 
obtain from 0 that 


z = (r 0 A) z -I- {Im ® B)w, (9) 

where T = diag{0, Iat-i} is a diagonal matrix including all 
eigenvalues of E in its diagonal. Next, let us partition U, the 
state and input vectors in 0 as follows. 


U= [Ui 


U 2 ] 




( 10 ) 


where Ui G rMx(m-7v+i)^ ^ ^ 

^n{M-N+l)^ 52 G tui G R"(J^-^+l)^ g 

R”(^“i). Then 0 is equivalent to 

zi = {Im-n+ 1 ® B)wi, 

h = {In-1 ® A)z2 + {In- 1 ® B)w2, 







3 


and zi = {Ui 0/„)z, Z 2 = (C/J ®In)z. Now, the initial MAS 
Q is transformed along the edges to a system composing of 
two independent subsystems. In the next section, we show that 
the consensus design for Q is equivalent to the stabilizing 
design for the second subsystem in (fTTl i. 


C. Design equivalence 

In this paper, we assume that all states of agents are 
measurable hence a distributed state feedback consensus con¬ 
troller for the MAS 0 is designed. However, if not all 
states of agents can be measured then we can employ a 
decentralized Luenberger-type observer as in cni to obtain a 
distributed output feedback consensus controller. Due to space 
limitation, we do not duplicate the design of such decentralized 
Luenberger-type observer in the current paper. 

Suppose that the stabilizing state feedback controller for the 
transformed edge dynamics (fTTl i has the form w = —Tz. Let 
r € be the diagonal matrix including all non¬ 
zero eigenvalues of L in its diagonal, and V € is an 

orthogonal matrix such that 


V'^LV 


0 0 
0 r 


( 12 ) 


Partitioning V into [14, V2] where 14 £ 14 G 

Then 

L14 = V 2 r ^ LV 2 = r, (13) 


since 14^14 = In-i- The following theorem shows the equiv¬ 
alence between the distributed stabilizing controller for the 
transformed edge dynamics (fTTl i and the distributed consensus 
controller for the initial MAS. 

Theorem 1: Suppose that Q is connected. Let U 2 be chosen 
as i5^14r“^/^, then the distributed stabilizing controller w = 
—Tz for the transformed edge dynamics with T = F ® K, 
K G R™^" and 


F = 


0 

0 


0 


(14) 


is equivalent to the consensus controller 


u = —n{L (g) K)x, 


(15) 


for the initial MAS. Furthermore, = 0 V f > 0. 

Proof: First, we show that the orthogonality of U is 
satisfied with U 2 chosen to be £'^14r“^/^. Indeed, C/J(7i = 
EUi = 0 since EUi = 0 due to a fact that 
EL = E. Furthermore, U^U 2 = EE^V 2 T-^^^ = 

p-l/2pp-l/2 ^ 

Next, multiplying to the left of (fTTl i with E"'" gives us 


L^E'^V2 = E^V2V, 

^ LeL;^i4r-i/2 = e'^V2V^/^, 

^ ELeE'^V2T-^/^ = p-^/'^V^EE'^V2P^/^, 

^U^LeU2 = P. (16) 


On the other hand, LUi = 0, which leads to LgC/i = 0 since 
LeL = Lg. Therefore, we obtain 


F = nU'^LeU. 


(17) 


This is equivalent to 

F{U'^E^) = nU^LeE'^ = fiU^E'^L. (18) 

Consequently, 

w = —Fz, 

^ [{U'^E'^) (g) /m]u = -{E g) K)[{U'^E'^) 0 In]x, 

^ [{U'^E^) 0 I^]u = -[{FU^E^) 0 K]x, 

[(W^E^) 0 Im]u = E"^L) 0 K]x. (19) 

Since E^ = [EUi,EU 2 f = [0,14r^/4'^, (O occurs 
if and only if [iP^^'^V^) 0 Im\u = — [pfF^/^l^^L) 0 K]x, 
which is equivalent to [{V 2 V 2 )®Im]u = — [/r(V 2 V 2 ^L) 0 A']x 
by multiplying to the left with (V 2 r“^/^) 0 Im- Note that 
14 = then 1414^ = Im — ViV^ = In — ^InIn- 

Hence, we obtain [(/at - ® Im]u = -[fi{[lN - 

^l 7 vl[^]L) 0 K]x = —fi{L 0 K)x, since iJfL = 0. This 
is equivalent to it = —IJ.{L 0 K)x -f (Iat 0 Im)uo for any 
uq G R"*. Since we are not interested in self-feedback inputs 
for agents, itg = 0 or equivalently u = —/i(L 0 K)x. 

On the other hand, we have Zi = z = [{Ui i7^)0/„]x = 
0 since EUi =0. ■ 

The interesting consequence of Theorem [T] is that the dis¬ 
tributed consensus design for 0 is equivalent to the synthesis 
of the distributed stabilizing controller W 2 = —{pP ® K)z 2 
for the transformed edge dynamics, which even has lower 
dimension. 

Remark 7.- If 0 is a spanning tree then M = N — 1 and 
hence L = In-i- Then we do not need the additional trans¬ 
formation 0. Therefore, all results here and in subsequent 
sections are derived with 'W 2 and Z 2 replaced by w and z, 
respectively. 

III. Full-Order Distributed Consensus Design 
A. Globally Optimal Stabilizing Design for Edge Dynamics 

Consider a performance index 

poo 

'^ = / {Z2Qz2+W2liW2)dt, ( 20 ) 

Jo 

where Q G R(N-i)nx(Ar-i)n^ Q A 0, G 

^{N-i)mx{N-i)m^ 7^ 0 are the weighting matrices. Sub¬ 

sequently, the stabilizing state feedback controller is 'W 2 = 
—E 2 Z 2 with 74 = IZ~^{Im-i ®B)'^V where V is the positive 
definite solution of the following global Riccati equation 

0 = F\Im—i 0 A] -|- [/at —1 0 F F Q 

-r{iN-i®B)n-\iN-i®Bfr. ( 21 ) 

It can be seen that (fTTfl is a distributed structured Riccati equa¬ 
tion and its solution V should also have a distributed structure 
such that the stabilizing feedback gain 74 has a specific form 
shown in Theorem [T] In the literature, there have been some 
studies on structured LQR designs for MASs, e.g. m-na, 
which assumed that the system matrices and the weighting 
matrices belong to some special classes, e.g. operator algebras 
or semigroups so that the solution of Riccati equation and the 
feedback gain also have that structure. Nevertheless, we do not 
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need those assumptions in our work. Therefore, our design is 
applicable for a broader class of systems. 

The following theorem reveals our design for distributed 
globally optimal stabilizing LQR controllers. 

Theorem 2: Let the weighting matrices be selected as fol¬ 
lows, 

Q = ®Qi + — In-i) ® Q2, 

where Qi e Qi h 0, i?i £ 

i?i 0, ^2 = PiBRiB^Pi, such that ,Im-i ® A) 
is detectable, where Pi is the solution of the following local 
Riccati equation 

PiA + A^Pi +Qi- PiBRiB'^Pi = 0. (23) 


Moreover, 


T > 


'^min (L)' 


(24) 


Then the unique solution of the Riccati equation (EB is 
V = /Ar_i(g)Pi. Furthermore, the optimal stabilizing feedback 
controller gain for the transformed edge dynamics is 


F2 = ^ir(^{RiB^Pi). (25) 


Proof: With the weighting matrices selected in (l22l i and p 
chosen in (l24l i. Q ^ 0 and 7^ 0. Note that In-i®A) 

is detectable if A) is detectable since it means (/jv-i® 

Qf , I]ii-i®A) is also detectable and adding the term (/rr — 
In-i) ® Qa ^ 0 does not break its detectability. Then the 
global Riccati equation (ISTTi can be rewritten as follows, 

0 =7^(/jv —1 0 A2j + (/^v —1 0 A)^P + In—i 0 Qi 

+ {fiT - In-i) 0 ^2 - P[pT 0 {BRiB^)]V. (26) 


Consequently, we can observe that (|26] | is satisfied if P = 
In 0 Pi, Q 2 = PiBRiB^Pi, and Pi is the solution of the 
local Riccati equation (|2B . Thus, V = In 0 Pi is, ^ solution 
of ( l26l l. On the other hand, the global Riccati equation (ISTT i 
has a unique solution, therefore V = In 0 Pi is indeed that 
unique solution. As a result, the optimal stabilizing feedback 
controller gain for the edge dynamics can be calculated as 
follows. 


F 2 =IZ ^{In-i 0 B)'^{In-I 0 Pi) = /iT 0 (PiP^Pi). 


Remark 2: The optimal performance index in this case is 

Af-l 

JnUn = 52(0)^P52(0) = ^ 'Z2 A^)^ Pl'^1 A^). (27) 

i=l 

which is independent of the information graph, where Z 2 (0) = 
[^ 2 , 1 ( 0 )^,..., -S 2 , 7 v-i( 0 )^]^ is the vector of initial conditions. 
Therefore, the cost for stabilization is not changed even if the 
communication structure in the MAS is varied. This is more 
advantageous than other studies, e.g. 0, where the minimum 
performance index depends on the communication graph. 

Remark 3: The condition for the coupling strength /r in 
Theorem |2] is similar to other existing results 0, 0, ifT^ . 
However, we recognize that this lower bound of the coupling 


strength is conservative, as will be completely removed in the 
next subsection and will also be demonstrated in the illustrative 
example. 


B. Locally Optimal Stabilizing Design for Edge Dynamics 

As seen from (fTTIi . the stabilizing design for the transformed 
edge dynamics is actually decentralized, which means that the 
stabilizing controller for the transformed edge dynamics can 
be designed locally, i.e., each edge is equipped with a local 
performance index as follows, with i = 1,..., N — 1, 



where Qi G Qi F 0, and Pi G Pi ^ 0 are the 

same as in Theorem |2] Then the controller gain for each edge 
dynamics is calculated by K = RiB^Pi, where Pi G 
Pi 0 is the solution of the Riccati equation (l23t . From 
LQR theory, it is known that the eigenvalues of A is shifted to 
the open left half complex plane by the LQR controller with 
its gain K. Therefore, by scaling with a parameter fiyt > 
0 V 1 < * < W — 1, the controller gain pyiRiB'^Pi still 
shifts the eigenvalues of A to the left though it could be more 
or less depending on whether p.yi > 1 or giyi < 1. Since 
we have assumed that all eigenvalues of A are on the closed 
left half complex plane, this means all eigenvalues of A — 
pyiBRiB"^Pi belong to the open left half complex plane for 
alH = 1,..., N—\. Thus, the transformed edge dynamics (fTTT i 
is stabilized by the locally optimal decentralized controller 
w = —Fz as in Theorem [T] and equivalently the initial MAS 
(0 reaches consensus with the controller (flsT l. for any p, > 0. 
This gives us the following theorem for the locally optimal 
distributed consensus design. 

Theorem 3: For any /i > 0, the transformed edge dynamics 
dnii is stabilized by the locally optimal decentralized con¬ 
troller w = —Fz as in Theorem [T] with the controller gain 
K = RiB^Pi and Pi is the solution of (|2B . Accordingly, 
the initial MAS 0 reaches consensus with the controller (ffSl l 
for any jj, > 0. 

The remarkability of Theorem [3 is that there is no other 
condition on the coupling strength p, rather than its positivity, 
which to the best of our knowledge has not been proposed 
hitherto. All existing results, e.g., 0,0,0, im required 
p to be greater than some certain bounds. Furthermore, those 
lower bounds of p are dependent on the graph Laplacian L, 
which is a global information and hence adaptive methods, e.g. 
ma, should be designed to overcome this drawback, or other 
methods should be developed to compute the lower bound of 
p in a distributed way. On the other hand, our result shows that 
no further techniques need to be implemented, which greatly 
simplifies both the theoretical study and the implementation 
of the consensus law. 

It is also interesting that the sum of all minimum local 
performance indexes (ESl l is ^ 2 ,i(0)^Pi 22 , 1 ( 0 ), which 

is the same with that of the global performance index in (l27l i 
since Qi and Pi in (|28] | are equal to those in (l22li and hence 
result in the same Pi, although the designs are different. 
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IV. Reduced-Order Distributed Consensus 
Controller Design 

A. General Approach 

Let (Ai,..., Ag) be eigenvalues of A that are unstable or 
needed to move further to the left of the imaginary axis and 
are the associated left-eigenvectors; 0 < q < n. We 


select the weighting matrix Qi under the form 

Qi = VQiV*, 


(29) 


where V = [i^i • • • i^q], Qi h 0. Denote 

A = diag{AJi=i,...,g, H = B^V. (30) 

We then show in the following theorem that this selection of 
Qi in addition to the designed consensus controller gives us 
a reduced-order consensus controller for the MAS Q. 

Theorem 4: By choosing the sub-weighting matrix Qi as 
in (|29] | such that C) A) is still detectable, the 

controller (flSl l designed in Theorem |2] or Theorem [3] has the 
order of at most q < n, resulting in a reduced-order distributed 
consensus controller for the given MAS ([^l. Furthermore, 
the eigenvalue set the closed-loop MAS matrix A can be 
characterized as follows, 

( U ^(= 7 )] U(®'(^)\'LAl,---,Ag}), (31) 

WGctIL) / 


where S-y is defined by 


= A - p-fH*RiHri. 


(32) 


Simultaneously, the eigenvalue set of the closed loop trans¬ 
formed edge dynamics can be defined by 


U 

< ■y^cr{L),-y^O 


lJ(a(A)\{Ai,...,Ag}). (33) 


ViA + AVi - ViTliVi -f Qi = 0, 


(34) 


Proof: Suppose that Vi € >- 0 is the unique 

solution of the following Riccati equation 


where 7^l = V*BRiB^V. Let Pi = VViVfPi >- 0, 
then substituting Qi and Pi back to the Riccati equation 
(l23]l, we can easily verify that Pi is a solution of (I23ll. 
Since ( |2^ has a unique positive semidefinite solution. Pi = 
VViV* is indeed that unique one. Since rank(PiP^Pi) = 
rank(PiP^VPi V*) < q, the controller designed in Theorem 
12] becomes a reduced-order distributed consensus controller 
whose order is at most q. 

Next, we will show that the eigenvalue spectrum of the 
closed-loop matrix A can be characterized by (l3Tl i. Consider 
any right eigenvector pj € C" of A corresponding to an 
eigenvalue Xj € cr(A)\{Ai,..., A,}, j = q + 1,... ,n. For 
i = 1,... ,N, denote Ci G the vector whose ith element 
is 1 while all other elements are zero. Then 

A(ei 0 rij) =ei 0 - p.{Lei) 0 (PPiP^PiPj), 

=Xj{ei 0 rjj) - p{Lei) 0 {BRiB'^V'PiV*r]j), 


due to a fact that V*rij = 0 since iy*r]j = 0 ^ I = 1,... ,q. 
Therefore, Xj is an eigenvalue of A with the associated right 
eigenvector Ci 0 rjj. This leads to 

C cr(A), (36) 

with the note that each eigenvalue Xj,j = q + 1,... ,n has 
the multiplicity N. On the other hand, 

{In 0 V*)A = {In 0 V*)[In ® A- p.L® (PPiP^Pi)] 

= In® {AV*) -pL® {V*BRiB^VViV*), 
= [In®A-pL®{H*RiHVi)]{In®V*). 

Let p* be any left eigenvector of L and 7 is the associated 
eigenvalue. Then multiplying both sides of above equation 
with p* 0 ^*, ^ S C^, we obtain 

{p* ®C){In®V*)K 

={p* ® C)[In 0 a - pP 0 {H*RiHVi)]{In 0 V*), 

= [p* 0 (r A) - P7P* 0 {CH*R^HVi)]{In 0 V*), 

=P* 0 (e*A - p7e* H*RiHVi){In®V*). (37) 

Now, suppose that is a left eigenvector of defined in 
(|3^ and the corresponding eigenvalue is a, it is deduced from 
(ETji that 

{p* ® n{iN 0 v*)A = p* 0 {e^j){iN 01^*), 

= a{p*®C){lNG)V*), 

= ap* 0 {CV*). (38) 

This means a is also an eigenvalue of A with the associated 
left eigenvector p* 0 Furthermore, a has the multi¬ 
plicity N. Accordingly, 

cr(S.y) C cr(A). (39) 

Thus, combining (l3^ and (l3^ gives us OTT i. 

The eigenvalues set of the closed-loop transformed edge 
dynamics can be proved in a similar manner. ■ 

B. First-order Distributed Consensus Controller 

Assuming that A has a single 0 eigenvalue. From the result 
of Theorem |4| we see that to ensure the consensus of the 
given MAS or the stability of the transformed edge dynamics, 
the reduced-order consensus controller should shift this 0 
eigenvalue of A to the left haft complex plane. Let G 
be the left eigenvector of A associated with the eigenvalue 0. 
Consequently, choose Qi = vqiu^ where qi > 0. Suppose 
that this selection of Qi still guarantees that Ijv-i0A) 

is detectable. The following result is immediately obtained as 
a corollary of Theorem |4| 

Corollary 1: The MAS (|3]l achieves consensus by a 1st- 
order distributed controller u = —p{L® K)x with any p > 0 
and 

. fnT „ „ 

(40) 


K = ^^RiB'^vv^, 


=Xj{ei ®r]j), 


(35) 


where ri = BRiB^v. The closed-loop eigenvalues of the 
MAS are 

U j U 

k7Go-(i) / 
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and of the transformed edge dynamics are 

IJ I IJ (ct(^)\{0}) . (42) 

7Gcr(L).75^0 / 

Furthermore, the consensus speed is determined by 

min{p,/giriAniin(i), Ai„in(-^)}- (43) 

The strong advantages of Corollary [T] are as follows. First, 
it gives us a superior simplicity for distributed consensus 
controller design since the controller’s order is just 1. Second, 
it clearly supports our relaxation on the coupling strength 
p > 0 in Theorem |3] by looking at the eigenvalue set (l42li 
of the transformed edge dynamics. And last, the consensus 
speed and the eigenvalue distribution of the closed-loop MAS 
are explicitly established, whereas they are implicit in the full- 
order distributed consensus controller design. In addition, the 
consensus speed is adjustable by varying p, qi, and ri. 

Remark 4: Although only the results of Ist-order dis¬ 
tributed consensus controllers are presented in this section, 
a similar design is readily obtained for 2nd-order distributed 
consensus controllers in the scenario that A has a pair of 
complex conjugate eigenvalues on the imaginary axis. Par¬ 
ticularly, Qi is chosen as in (|29] | with V composes of two left 
eigenvectors associated with the eigenvalues on the imaginary 
axis of A. 

C. Illustrative example 

Consider a drying section of a paper converting machine 
including of 9 rolls Eol where the control target is to syn¬ 
chronize the angles of the rolls in order to make the paper run 
smoothly without tearing. Assuming here that each roll just 
communicates with the rolls in front of and behind it. The 
state matrices of the fcth roll, fe = 1,..., 9, is 


'0 

1 

0 


'O' 



0 

-0.01 

0.2 

,B = 

0 

,C = 

0 

0 

0 

-125 


20 


0 


It can be seen that A has a single 0 eigenvalue while the other 
two eigenvalues belong to C_. Therefore, we apply the Ist- 
order consensus controller design in Corollary [T] with i?i = 
100 and qi = l. Figure [T] displays the consensus of the rolls’ 
outputs as p = 7 for globally optimal design, and p = 0.01 for 
locally optimal design. We can observe that even the coupling 
strength p is very small but positive, the consensus is still 
achieved. This confirms our result on the relaxation of the 
conservative condition for the coupling strength p,. 

V. Conclusion 

This paper has proposed a systematic approach to design 
distributed consensus controllers for general linear MASs with 
the following remarkable features. First, a transformation has 
been introduced, based on which the distributed consensus 
design is shown to be equivalent to a distributed stabilizing 
design for the transformed system called edge dynamics. Then 
the conservative bound for the coupling strength is completely 
removed by the locally optimal synthesis, resulting in a fully 




Fig. 1. Consensus in the drying section of a paper converting machine by a 
Ist-order consensus controller with different values of the coupling strength. 

distributed consensus controller. Second, a method has been 
presented to reduce the consensus controller’s order which can 
be 1 if the agents have a single pole at the origin. This gives 
us not only a superior simplicity for the distributed consensus 
synthesis, but also an explicit description of consensus speed 
and eigenvalue distribution of the MAS. 
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